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Abstract. IceCube is a cubic -kilometer neutrino telescope under construction at the geographic 
South Pole. Once completed it will comprise 4800 optical sensors deployed on 80 vertical strings at 
depths in the ice between 1450 and 2450 meters. Part of the array is already operational and data was 
recorded in the configurations with 9 (year 2006/2007), 22 (year 2007/2008) and 40-strings (year 
2008/2009) respectively. Here we report preliminary results on the search for point-like neutrino 
sources using data collected with the first 22 strings (IC-22). 

Keywords: Cosmic neutrinos, high-energy gamma-rays 

PACS: 95.55.Vj, 95.55.Ka, 95.75.Mn, 95.75.Pq, 95.80.+p, 95.85.Ry 

RATIONALE OF HIGH ENERGY NEUTRINO ASTRONOMY 

High energy neutrino telescopes attempt to resolve the origin of cosmic rays by de- 
tecting neutrinos. Whatever the sources of cosmic rays are, they are expected to emit 
high energy neutrinos as well. While hard to detect, neutrinos have the advantage of 
representing unique fingerprints of hadron interactions and, therefore, of the sources of 
cosmic rays. Moreover, their directions are not scrambled by magnetic fields. 

The candidates targeted by high-energy (E> 1 TeV) neutrino telescopes are astrophys- 
ical objects showing clear evidence of non thermal emission and providing hints for 
cosmic ray acceleration. A class of most favored candidates for cosmic ray acceleration 
up to the "Knee" are supernova remnants (SNRs), while only two types of objects are 
considered capable of accelerating protons at higher energies (up to 10 20 eV): gamma- 
ray-bursts (GRBs) and active-galactic-nuclei (AGNs). These offer the brightest events in 
the universe but mostly lie at cosmological distances, hard to observe with both protons 
and photons. While high energy cosmic rays interact with the cosmic microwave back- 
ground and photons with the extragalactic background light (currently the maximum 
distance tested with observed high energy gamma-rays is z=0.536 [1]), neutrinos would 
allow to extend the observable horizon of the universe. 

High energy gamma-rays are successfully detected for two decades from a large 
variety of sources (galactic and extragalactic) using satellites and air shower detectors. 
Moreover, the energy spectra can now be investigated in detail by new generation 
Cherenkov telescopes, like H.E.S.S., MAGIC and VERITAS. However, the observation 
of high energy gamma-rays alone does not provide unambiguous proof of hadron (i.e. 
cosmic rays) acceleration, since in most cases, the observed fluxes and spectra can be 
ascribed to high energy electrons and positrons which, once accelerated, give rise to 



gamma-rays either via inverse Compton scattering or bremsstrahlung. 

In parallel Ultra high energy cosmic rays astronomy, is hunting for the sources of 
cosmic rays at energies larger than about 10 19 eV, where the gyroradius of a proton 
in the galactic magnetic field is larger than the thickness of the galactic disk and the 
trajectory deflections are reduced. However, the cosmic ray flux at these energies is so 
small that the statistics so far collected does not allow us to pinpoint their origin. The 
Auger detector is designed to collect sufficient statistics to address this question and its 
recent results provided evidence of anisotropy in the arrival directions of the highest 
energy cosmic rays (at energies larger than about 60 EeV [0]). However, the possibility 
of identifying individual sources is debated (see for example [SIH]). 

Several cases of high energy gamma-rays from SNRs have been recently discovered 
(e.g. $$\). A remarkable example is RX J1713. 7-3946, for which observational data favor 
the hypothesis of cosmic ray (proton) acceleration (e.g. the hadronic model in [60 as 
opposed to the leptonic model in [70). X-ray measurements by Chandra, combined with 
detection of gamma-rays of energies up to at least 100 TeV, provide indeed a strong 
indication for acceleration of protons to energies of 1 PeV and beyond [8]. This source 
and similar ones may give measurable fluxes of high-energy neutrinos as well, if gamma- 
rays are generated by the decay of neutral pions (see for example |9l LlOQ for neutrino 
fluxes from SNRs and [l 1Q from the Cygnus region). Binary systems are also promising 
sources of neutrinos (e.g. lfl2ll for LS 5039), but various scenarios are predicted for 
the observed high-energy gamma-ray emission, each with different implications for 
the associated neutrino signal (e.g. 1130). Other objects also provide a good target for 
neutrino observations, whenever the observed gamma-rays have energies exceeding tens 
of TeV, like MGRO J 1908+06, observed by H.E.S.S. up to about 40 TeV [Q and by 
Milagro up to about 100 TeV [15]. The associated neutrino fluxes may be detectable by 
IceCube (e.g. 

Among extragalactic sources, an intriguing case is the blazar 1ES 1959+650, which 
manifested in 2002 an exceptional enhancement in its high energy gamma-ray emission 
which was not accompanied by an increase in the X-ray fingerprint of the entity expected 
in most leptonic scenarios [17Q . This so called "orphan flare" is considered an hadronic 
indicator. AMANDA reported high energy neutrinos in possible temporal coincidence 
with this period of exceptional activity, but this result is not statistically significant ||18|1 . 
Hadronic interpretations of this flaring episode also did not yield clear conclusions on the 
genuine nature of the detected neutrinos as a possible cosmic signal [T^, 20]. In order to 
reach the sensitivity required by most astrophysical scenarios to detect cosmic neutrinos, 
at least cubic-kilometer-scale arrays are required 112 10 . For models predicting high energy 



neutrinos in AGNs see for example [220 (and references therein), and [|23L I24Q for 
neutrinos from GRBs. 

Large data samples from AMANDA-II B25Q have been analyzed, showing no evidence 
of point-like signals from the northern sky in either time-integrated or time-variable 
searches [IflEI.HZ]. In [26], the sensitivity to tau neutrinos was estimated for the first 
time. An improved analysis based on a maximum likelihood method was later developed 
and an analysis of the complete AMANDA data set (about 3.8 years of effective lifetime) 
has been presented, yielding the most stringent flux upper limits to point sources of 
neutrinos 11270 . Results based on IceCube data (9 strings) have also been shown. No 
indication of point sources has been found yet [270. 



Neutrinos in association with GRBs are also being searched for in dedicated analyses. 
Besides promising sources of neutrinos, they offer the advantage of higher signal-to- 
noise ratio compared to steady point-sources due to the very small burst time window 
in which the neutrino signal is expected. Since IceCube has a large field-of-view, a very 
large portion of the sky is simultaneously observable and data can be scanned to look 
for correlations with GRBs catalogs. The best limits so far obtained are from AMANDA 
data. Optimistic predictions are ruled out ll28ll . 

New search strategies are recently being developed, aiming at increasing the chance 
to discover comic neutrinos by on-line searches for correlations with established signals 
(e.g. flares in high-energy gamma-rays) triggered by neutrino observations (neutrino- 
triggered target-of-opportunity). For sources which manifest large time variations in the 
emitted radiation, the signal-to-noise ratio can be increased by limiting the neutrino ex- 
posures to most favorable periods. The chance of discovery can then be enhanced (the so- 
called "multi-messenger approach") by ensuring a good coverage of simultaneous data 
at a monitoring waveband (e.g. gamma-rays). The first realization of such an approach 
led to two months of follow-up observations of AMANDA triggers by MAGIC 02911 . An 
extension of this program to IceCube and also to optical follow-up observations is in 
progress. 

A search for a cosmic neutrino signal is finally performed by looking for distortions in 
the background properties at high energies (above TeV), where the flux of atmospheric 
neutrinos is much reduced. Also in this case only upper limits are so far derived (cfr. 
||30|l for limits on diffuse fluxes from unresolved sources across the northern sky based 
on AMANDA data in the energy range of 16 TeV to 2.5 PeV, 113 ill for limits across the 
horizon in the energy range of 0.2 PeV to 1 EeV and for all-sky limits from cascade-like 
events Il27ll ). Results from IceCube (nine strings) have been reported in 

Besides hunting for astrophysical neutrino sources, IceCube is dedicated to particle 
physics topics including the indirect search for Dark Matter, the study of Lorentz Invari- 
ance, the study of neutrino oscillations and the search for magnetic monopoles. 

Here we report the latest results on the search for neutrinos from point-like sources, 
as obtained from the analysis of data collected with IceCube (22-strings configuration), 
with an effective lifetime of 276 days. 



THE ICECUBE TELESCOPE & OPERATION PRINCIPLE 

IceCube is a cubic-kilometer neutrino telescope under construction at the South Pole. 
Once completed in 2011, it will comprise 4800 Digital Optical Modules (DOM) along 
80 vertical strings deployed in the ice at depths from 1450 to 2450 meters. The first nine 
in-ice strings (IC-9) recorded data in 2006-2007 ||3~2|l and the first 22 in-ice strings (IC- 
22) from May 2007 to April 2008. Forty strings (half array) are operational since May 
2008. 

The neutrino telescope is complemented with an air-shower array (IceTop) which 
covers a surface of one square kilometer and consists of 80 detector stations, each 
equipped with two ice Cherenkov tanks. IceTop serves as a veto to cosmic ray showers 
and as a calibration array, and covers in addition science topics on cosmic rays physics. 
In 2006, the particle spectrum of the solar flare of December 13 was measured in the 



energy range between 0.6 and 7.6 GeV rf33ll . 

IceCube is optimized to detect neutrinos of all flavors in the energy range from TeV to 
EeV. Signal waveforms from individual photomultipliers are digitized in each DOM and 
transmitted to the surface by twisted-pair cables. Data is processed by several hardware 
and software triggers, responsive to different event topologies. A series of reconstruction 
algorithms and event filters are applied in situ to calibrated data to reduce the trigger rate 
(about 650 Hz for the SMT with 22 IceCube strings) to a total data volume of about 30 
GB/day, which was then transferred north via satellite. Events are then reconstructed 
using more CPU intensive maximum-likelihood techniques, in which events are fitted to 
templates representing different neutrino interaction modes. For example, and v z CC 
interactions can lead to muon tracks (in case of v T arising from the muonic decay mode 
of the tau lepton), which, at TeV energies, provide a good handle on the primary neutrino 
direction. On the other hand, v e and v T CC, as well as all-flavor NC interactions, can lead 
to compact energy deposits (cascades, in case of v T CC interactions arising from both 
the hadronic and electronic decay modes of the tau lepton), which provide a good handle 
on the primary neutrino energy. 

RESULTS ON POINT-LIKE SOURCES OF NEUTRINOS WITH 22 

ICECUBE STRINGS 

Searches for astrophysical sources of neutrinos have to cope with background events 
from the interaction of cosmic rays in the Earth's atmosphere. Decays of secondary 
mesons induce downward going muons and a more isotropic flux of neutrinos. Rates 
measured with 22 IceCube strings are O(10 9 ) events/year from downward-going atmo- 
spheric muons and 0(1O 3 ) events/year from atmospheric neutrinos. In view of point- 
source searches, neutrino candidates are typically selected starting with angular cuts, 
i.e. rejecting downward-going muon tracks, since only neutrinos can travel through the 
Earth. This limits the sensitivity to the northern sky. An effort is on-going to overcome 
this limit and extend both the field of view and the energy window of sensitivity (i.e. 
multi-PeV) & 

Three independent analyses developed for point source searches with 22 IceCube 
strings in the energy region from TeV to PeV are presented here: a binned search, based 
on an event counting approach using circular angular search bins, and two un-binned 
searches, one applied to IceCube data and based on a maximum likelihood approach 
which models the signal expectation based on the detector point spread function and a 
simple energy estimator, and another optimized for soft spectra, applied to combined 
IceCube and AMANDA events. 

The event selection is based on similar parameters in all three analyses. First, a 
filter based on event direction and quality parameters is used to reduce atmospheric 
muons by about three orders of magnitude. Yet an irreducible background remains 
from upward-going muons induced by atmospheric neutrinos together with a small 
fraction of mis-reconstructed downward-going muons plus possible signal events. A 
signal from cosmic neutrinos would manifest itself as an excess of events over the 
background at a specific location in the sky. Further track quality cuts are therefore 
applied, designed to achieve a good angular resolution, including angular uncertainty 



estimates, the number of un-scattered photon hits with a small time residual w.r.t. the 
Cherenkov cone, a likelihood ratio comparing a prior of an atmospheric muon track 
(downward-going) to the reconstructed one (upward-going), the reduced likelihood of 
the directional reconstruction, etc. To avoid biases in the event selection, the final cuts are 
optimized in a blind approach. In the search for point sources, where the event direction 
is used to look for a signal, this is accomplished by fixing cuts on a sample of events 
with randomized right ascension. In a given declination band (of size comparable to the 
angular resolution), a negligible point source neutrino contribution compared to the total 
background is expected, allowing an estimate of the background from scrambled data. 
Details of each analysis and respective results are given below. 

The IC22 binned analysis distinguishes a localized excess of signal from a uniform 
background using circular angular search bins. The sensitivity (average upper limit 
in case of no signal), in a given declination band, can be obtained from Poisson 
statistics and it is calculated by comparing the number of simulated signal events 
for a given flux reconstructed inside the search bin, n s j g , to the average number 
of background events as measured from scrambled data ni ?g , using Il26ll : 

® lim = ® Q - m {n bg )/n sig {<$>°) (1) 
where /I90 is the event upper limit at 90% CL obtained in the absence of signal 



following the Feldman&Cousins ordering principle |34|]. Cuts on track quality pa- 
rameters were optimized considering the detection and reconstruction efficiencies 
as a function of declination. The search bin radius was optimized as an additional 
free parameter (mean value = 2.1°). Final cuts were chosen to deliver an optimum 
sensitivity for a compromise between a spectral index 7 = 2 and 7=3 for the neu- 
trino energy spectra. This analysis uses a different event directional reconstruction 
then the other two analyses, that delivers better angular resolution at very high en- 
ergies. However, it uses reduced event information (events are equally weighted, 
disregarding the individual directional error and energy) and it is therefore ex- 
pected to achieve a lower sensitivity. 2956 events have been selected. The atmo- 
spheric neutrino content of the samples is greater than 90%. The rest comes from 
mis-reconstructed down-going muons, mostly multiple events from independent 
cosmic ray showers. From simulation, a sky- and energy-averaged median angular 
resolution of 1.3° is estimated for signal neutrinos with spectral index 7=2. 
The IC22 un-binned analysis constructs a likelihood function which depends on the 
signal probability density function (pdf), obtained from MC simulation of the 
event energy and angular distribution, the background pdf (obtained from data as a 



function of the event position x/), and total number of data events n tot [351] 



i—ntot n n 
L = H —Si{xt,x s , Ei, 7) + -Zt-Biixi, Ei). (2) 

i= l n tot n tot 



The signal pdf 5, is characterized by each event's direction, x,-, angular uncertainty 
around x s , energy estimation, Ei (i.e. the number of DOMs with at least one hit), 
and the assumed energy spectrum. The number of signal events, n s , is found by 
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FIGURE 1. Left: Effective area of IceCube with 22 strings versus energy and declination. The decrease 
at high energies for vertical tracks is due to absorption of neutrinos in the Earth. Right: Sensitivity 
(average upper limit in case of no signal, at 90% CL) and discovery potential (minimum flux delivering a 
5(7 detection with 50% probability) for the binned (blue) and un-binned (red) analysis. The latter shows 
better performance, due to event weighting including the angular reconstruction error and a simple energy 
estimator. 



maximizing the likelihood ratio of the background plus signal hypothesis against 
the background-only case. 51 14 events have been selected. From simulation, a sky- 
and energy-averaged median angular resolution of 1.5° is estimated for signal neu- 
trinos with spectral index y = 2. Based on signal simulation with a spectral index 
7=2 this analysis provided both the best discovery and limit setting capability and 
was therefore a-priori chosen to deliver the final significances and upper limits. The 
enhancement gained by using the un-binned method (without the systematic error) 
is 35%. It also shows on average a factor of 2 improvement over that of the total 
statistics collected by AMANDA-II rf27ll and provides the best limits to date. 

The final effective area of the binned analysis is presented in Fig. [Hand is similar for the 
two analyses (binned and un-binned). Sensitivities and discovery potentials, as functions 
of declination, are also compared in Fig. [TJ The sky-average sensitivity at 90% CL to a 
generic flux (with spectral index y = 2) is <& V)1 - E 2 <1.3 (2.0) -10~ n TeV cm _2 s _1 for 
the un-binned (binned) method respectively. The statistical overlap of the event samples 
tested by the two analyses is 72%. 

A pre-defined list of 28 sources (galactic and extragalactic) gave no significant excess 
over the expected background. The biggest deviation was found for 1ES 1959+650 (Crab 
Nebula) with a post-trial p-value of 66% (91%) by the un-binned (binned) analysis. 
Details are given in Tab. [TJ 

The results from the all-sky search for both analyses are shown in Fig. [2] and Fig. [3] 
The maximum positive deviation from the background hypothesis was found at right as- 
cension a =153° (183.5°), declination 5 = 11° (47°), with a post-trial p-value of 0.67% 
(67.6%) for the un-binned (binned) analysis respectively. The post-trial p-value for the 
highest excess is 1.3%, taking into account a trial factor of two arising from the mul- 
tiple tests performed (i.e. pre-selected catalog and northern sky un-biased search). The 



TABLE 1. Pre-trial p- values from the 
binned and un-binned analysis of a pre- 
defined catalog of 28 sources (preliminary). 
Sources were selected based on evidence for 
non-thermal emission. In bold are marked the 
respective highest excesses. 
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spot found by the un-binned analysis (a-priori chosen to deliver the final significances) 
was further investigated (e.g. by a time-dependent analysis) resulting in no significant 
event cluster in time (no events are closer together than 10 days). It will be tested with 
better sensitivity using data with 40 strings. The difference in the p-value found by the 
two analyses arises mainly from the energy parameter implemented in the likelihood 
function. 

The IC22+AMANDA low-energy optimized analysis is designed for galactic 
sources, which are expected to exhibit softer spectra than a simple power-law 



FIGURE 2. Pre-trial significance sky map from the binned analysis (preliminary). 




FIGURE 3. Pre-trial significance sky map obtained with the un-binned analysis (preliminary). 



with 7= 2. A likelihood function, similar to the one given above, is implemented: 

f=»fO! n n 

L = EI —Si + (1-—)B,. (3) 

i= l n tot n tot 

with constant background pdf B[ = 1 /CI and a signal pdf 5,- = exp( — y 2 / (2a; 2 ) ) / 2n ■ 
<7j 2 ), with \j/ being the angle between source sky position and each event's recon- 
structed direction. All events in a ±5° band around the galactic plane (covering a 
solid angle CI) are included, with a lifetime of 276 days. 8727 events have been 
selected. From simulation, a sky- and energy-averaged median angular resolution 
of 1.5° is estimated for signal neutrinos with spectral index 7= 2. 
Four pre-selected sources (Crab Nebula, Cas A, SS433, LSI+61 303) were also 
tested with this method. 

This analysis also found no significant deviation from the background hypothesis, with 
the highest excess in the galactic scan found at longitude / = 75.9°, latitude b = 2.7° 
(post-trial p-value 95%) and the highest excess among the four pre-selected sources 
found in the direction of the Crab Nebula (post-trial p-value 37%). 



We conclude that no significant deviation from the background hypothesis was found 
in the IceCube data collected with 22 strings by the three analyses presented here. All 
results shown are preliminary. Flux upper limits, including the systematic uncertainties, 
are being calculated. 

SUMMARY & OUTLOOK 

The IceCube neutrino telescope is now more than half-way completed with 40 strings 
in operation in 2008. We have shown preliminary results on the search for point- like 
signals in the northern sky, using data collected from May 2007 to April 2008 with 22 
strings. No significant excess of events was found by three independent analyses. An 
additional work, focused on extending the field of view to the southern sky for multi- 
PeV neutrinos, will be presented elsewhere soon. Analysis of data in the 40 strings 
configuration is on-going. 

Data collected with 22 strings is also being searched for an overall excess due to un- 
resolved sources (diffuse flux). This work is expected to exceed the current best limits 
obtained with AMANDA-II (4 years of data), due to improved energy estimation and 
the larger acceptance above 100 TeV. 

The full IceCube neutrino telescope is expected to achieve a track angular resolution 
of less than one degree and an order of magnitude improvement over the sensitivity to 
point-like neutrino sources obtained with the entire data set collected with AMANDA- 
II (7 years of data), within three years of operation. Based on current estimations, the 
cubic-kilometer scale neutrino telescope has good chances to identify fluxes from steady 
point sources of high energy neutrinos. One year of data with 80 strings is expected to 
achieve a sensitivity to diffuse fluxes of neutrinos currently predicted from unresolved 
AGNs. 

Detection of high energy neutrinos with the next generation neutrino telescopes (Ice- 
Cube, KM3NeT) will probe the most powerful accelerators and the physical mechanism 
leading to their non-thermal emission. Increased discovery potential and a more compre- 
hensive phenomenological understanding is expected by combining neutrino data with 
established non-thermal astrophysical signals (e.g. high energy gamma-rays in the so- 
called "multi-messenger approach"). 
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